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Abstract. In this paper, we derive the Ewald method for inverse power-law 
interactions in quasi-two dimensional systems. The derivation is done by using two 
different analytical methods. The first uses the Parry's limit, that considers the 
Ewald methods for quasi-two dimensional systems as a limit of the Ewald methods 
for tridimensional systems, the second uses Poisson-Jacobi identities for lattice sums. 
Taking into account the equivalence of both derivations, we obtain a new analytical 
Fourier transform intregral involving incomplete gamma function. Energies of the 
generalized restrictive primitive model of electrolytes (?/-RPM) and of the generalized 
one component plasma model (77-OCP) are given for the tridimensional, quasi-two 
dimensional and monolayers systems. Few numerical results, using Monte-Carlo 
simulations, for 77-RPM and 77-OCP monolayers systems are reported. 
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1. Introduction 

A quite general class of long range potentials is the inverse power law potentials. 
Generally, these potentials are defined as pseudo-potentials or as effective potentials ; the 
lattice sums with these potentials are used in solids state physics for the computation of 
structural integrals [IJ|2l[3llIJ[5j|6j[7J|8] ; for instance, in the computation of the energy 
needed to permit the formation of an atomic vacancy in metals [6] or also to study 
the effect of a piezoelectric medium on structural properties of electronic bilayers in 
heterogeneous junctions AlGaAs-GaAs [9]. The dependence of the liquid-gas transition 
for inverse power law interactions with 77 = 3 + a has also been studied in reference [10] . 
Although, the analytical form of pseudo-potentials or effective potentials used are more 
complicated than a simple inverse power law, but rather like f(r)/r v , in the following, 
we will restrict ourselves to f(r) = 1. 

For a system with periodic boundary conditions in the three dimensions of the space, 
lattice sums of pseudo-potentials are 

= £ — l jwz } <u 

T (d) I 7 + I 

with Ln a symbolic notation for the d dimensional periodic images of the basic cell that 
contains the system ; the basic cell is noted Lq (or Lq for tridimensional periodicities, 
or Sq for two dimensional periodicities). 

For rj — 1, this class of potential corresponds to the Coulomb interactions which lattice 
summations have been derived in various ways including Ewald methods [HI [12], Lekner 
[TBI [HI [15] or plane- wise summations [T6l \T7\ . 

The case rj = 3 is also of a large interest in liquids and condensed matter physics, 
especially when dipolar interactions between particles are present in the systems [TT1[T8] . 
The limit 77—7-00 corresponds to the hard sphere model ; indeed, one has 

lim — = +00 for r < 1, and lim — = for r > 1. (2) 

rj— >-oo rj— >oo 

The limit r; — ?- corresponds to an interaction potential that do not depend on the 
distance between the particles (lim^o r v = 1) ; thus, from the point of view of computer 
simulations, the limit 77 — > corresponds to the ideal monoatomic gas. Logarithmic 
interactions or coulomb interaction in two dimensional systems cannot be obtained in 
the limit 77 — > of inverse power law interactions of Eq.flTJ). However, it may be obtained 
with the help of the relation 



lim 

jj->-0 



— lnr (3) 



The Ewald method for logarithmic interactions has been derived by Perram and de 
Leeuw [19]. One has to take into account Eq.([3]) to recover their results from the Ewald 
sums for inverse power-law interactions developed in this paper. 

Inverse power law interactions are also used in long-range dispersion force in Lennard- 
Jones fluids (77 = 6) ; some significant corrections are obtained by using Ewald 



Ewald sums for inverse power-law interactions 



3 



summations in Lennard- Jones fluid rather than using a truncation of the potential 
[201 EH [22]. 

The lattice sum (JTJ for arbitrary n has been obtained by Misra, Born and Bradburn 
[H [2] for tridimensional systems. The main purpose of the present paper is to compute 
this lattice sum for quasi-two dimensional systems (systems with spatial periodicities in 
only two directions of the space). 

There are mainly two ways to obtain Ewald methods for quasi-two dimensional systems. 
First, one may obtain this formulation by taking the so-called Parry's limit [23]; it 
consists in assuming that one of the spatial periodicity (say L z ) is extremely large 
compared to the other, the Ewald summations for tridimensional systems in the limit 
L z — > oo are Ewald summations for quasi-two dimensional systems. This derivation has 
been applied successfully to Coulomb, Yukawa and dipolar interactions [T8l [24] . The 
second method consists in the application of Poisson-Jacobi identities (discrete Fourier 
transforms) to the two dimensional periodicities of the quasi-two dimensional systems. 
Both derivations are fully equivalent and they provide exactly the same results, this was 
shown extensively for Coulomb, Yukawa and dipolar interactions [24" 1 [25l [T8] . 
Another method is based on a decomposition of the charge distributions with the help 
of different screening functions [261 EH [28j [291 USB EI]- This method has been applied to 
coulomb interactions [261 12H [28], [29] , to Yukawa potentials [30] and inverse power-law 
interactions [31]. The screening of charges is frequently taken as a gaussian function, 
but this choice is not restrictive [261 [2H [29] . 
The Poisson-Jacobi identity for lattice sums reads as 



where Vd is the volume of the basic cell, IZd the reciprocal lattice associated with the 



In this paper, we will derive Ewald summations for inverse power-law interactions in 
quasi-two dimensional systems by using both methods ; then, using the equivalence 
between both derivations, we will obtain a new analytical relation for a Fourier transform 
involving the incomplete gamma function. This Fourier transform is expressed in term 
of the generalized incomplete gamma function or incomplete Bessel function [321 [33] . 
The paper is organised as follows. In the next section we give a short derivation of 
the Ewald method for the lattice summations with inverse power-law interactions in 
three dimensional systems. Then, using the result of section 2, in section 3 we derive 
the Ewald method for quasi-two dimensional systems by taking the Parry's limit and 
also by using the Poisson-Jacobi identity ; the equivalence between both derivations 
of the Ewald method for quasi-two dimensional systems then allows us to obtain a 
new Fourier transform for incomplete gamma function. In section 4, energies of the 





lattice made of the periodic images of the system and / the Fourier transform of / 
defined as 




(5) 
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generalized restrictive primitive model of electrolytes (77-RPM) and of the generalized 
one component plasma model (77-OCP) are given for the tridimensional, quasi-two 
dimensional and monolayers systems. Several numerical results, using Monte-Carlo 
simulations for ?7-RPM and 77-OCP monolayers systems are reported. 



2. Tridimensional systems. 

To apply the Ewald method to inverse power-law potentials, we begin as it is done for 
Coulomb interaction, by using the relation 

1 i r dt 



Ln\'-T®h tM) *M-\r + L n \ t) (6) 

then we split the lattice sum into two summations using a convergence parameter a. 
One of the summations is on the periodic images of the system and the other one is 
transformed by using Poisson-Jacobi identity, both being rapidly convergent. We have 



Ln, v 1 ' 



V 1 f°° dt c -\r+L n \ 2 t , 1 r dt y> \ r +L n \ 



(7) 



For the second contribution, we use the Poisson-Jacobi relation written as 
we found 



Ln ken 



i„(r) = (V, "I I r + In |) + ^ $l 3) (r/, a; fc) exp (jfc.r) 



(9) 



7T" 



3 / 2 r 2 ^ 



r (2) v y t( 5 -^/ 2 

where the last term is the contribution for k = (it diverges for < 77 < 3 and it is finite 
for 77 > 3). For Coulomb interaction, the contribution for k = is closely related to 
the macroscopic boundary condition : a supplemental boundary condition imposed very 
far from the basic cell [121 El- Without further specification, the k = contribution is 
diverging ; this can be cancelled in coulomb interaction by assuming electroneutrality of 
the system in the basic cell. Similarly, for inverse power-law potential this contribution 
is diverging if 77 < d (d = 3, in this section - in section 4, we show how electroneutrality 
of the system in the basic cell suppresses this divergence). The way by which this 
IR-divergence is cancelled for inverse power-law potentials depends on the particular 
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Power Real Space 

Contributions $ 



(3) 



Reciprocal Space 
Contributions <£>, 3 ^ 



Contributions 
for k = 



r) — 4 

77 = 3 

77 = 2 

t; = 1 



'i + a V) 



r / \ i 2ar -a 2 r 2 

crlc(ar) H j=e 

V 7r 



crfc(ar) 



7T 



— fc erfc ( — 
2a 



~ Ei(-fc 2 /4« 2 ) 



2?r 2 erfc(fc/2a) 



2tt 3 / 2 

( 

V 



Am , / a \ 



2^V2 
V 



1 1 

e a 



1 1 



Table 1. Analytical formulas for the real and reciprocal space contributions in 
tridimensional Ewald summations for inverse power law interactions 1/V for 

(3) (3) 

some integer values of rj. & R and $L contributions are given by For 
n < 3, the contribution for fc = has a diverging behaviour which asymptotic 
expansions are given by using the prescription e defined in equation ([TO]). 



physical situation and on systems one studies. This IR-divergence is outlined below by 
the introduction of the prescription e as 

' alt .. r 2 alt 



£(5-77)72 



lim 



£(5-77)72 



(10) 



that gives the asymptotic expansion of the last integral in Eq. 
Both functions and $/<; are given by 



r g) 



$f foa; *) 



7T 3 / 2 /_4 



(3-g) F / 3-?7 fc 2 \ 
2 1 ^~2~' 4^ J 



(11) 



v rg) 

with T(a,z) the complementary incomplete gamma function. These results agree with 
the derivations done in refs. [20], [21] ; Table [T] gives some analytical forms of real and 
reciprocal contributions for few integer values of r\. 



3. Quasi-two dimensional systems. 

To compute surface properties based on effective or pseudo-potentials, Ewald 
summations for quasi-two dimensional systems are of interest [9]. Quasi-two dimensional 
systems are heterogeneous systems with some anisotropics in their spatial extensions, 
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their numerical studies are done with partial boundary conditions : periodic boundary 
conditions are taken in directions with large spatial extensions and other boundary 
conditions are taken in directions with smaller extensions. For these systems, the lattice 
sums of inverse power law interactions are given by 

( r ) = y (12) 

^ \ r + Sn\ r ' 

(2) 

with Sn (= Ln) the symbolic notation for periodic images of the system in the 2D 
geometry. 

As stated in the introduction, there are mainly two analytical ways to derive the Ewald 
method for a quasi-two dimensional system with a given interaction. First, if one already 
knows the analytical form of the Ewald summations for the corresponding full periodic 
tridimensional system, then one may make particular one direction of the space (Oz 
axis for instance) along which no periodic image is taken ; this is achieved by taking the 
limit L z — > oo with L z the periodicity in the z direction of the tridimensional system. 
This derivation was obtained by Parry for Coulomb interactions [23] and also, more 
recently, for Yukawa interaction [24] ; it can also be done for dipolar interactions. For 
systems in the Parry's limit, we use the notations 

V = L Z A ; r = s + ze z ; L n = S n + n z L z e z and k = G + k z e z (13) 

with the component k z of the wave vectors in the Oz-direction written as k z = 2irn z /L z 
with n z integer and the two dimensional wave vectors G belonging to the reciprocal 
lattice associated with the two dimensional lattice defined by the periodic boundary 
conditions. 

When the limit L z — > oo is taken, the reciprocal space contribution has to be considered 
with caution. According to notations of equation (1131) . one has to separate the 
summations on the wave vectors as 

• = '* : W if G ^ °> then kz = is allowed ; 

k^o G (14) 
(ii) if G = 0, then k z ^ 



+oo 

k=io G^o 



With these notations, the summations on the reciprocal vectors are split as 

1 



where the contribution Tq_^ corresponds to the case (ii) of Eq. (fT4"]) . it is given by 



(15) 
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and it includes a non trivial contribution that depends on z and also, if 77 < 2, it has an 
IR-divergence. 

The summations over periodic image in the real space are also split as 



L n S n 



r(f,« 2 


| (s + S n ) + z&z 


I 2 ) 


1 


(s + Sn) + ze z | 





;i7) 





| (s + 5 n ) + + n z L z )e z 


1 2 ) 


r(l)l 


(s + S n ) + (z + n s L z )e z | 


V 



In the Parry's limit (L z — > 00), only the first contribution in right hand side of Eq.( TT7|) 
survives, then Eq. ffl2l can be cast into the form 



(v) 



'a,z) (11 



</>„(r) = ® { r (v,a;\s + S n + ze z \)+J2 $ g % «, Z] G)e> G - s + T ( 
S n G+q 

where $^(77, a; r) is given by Eq.ffTB and $[? 2) (r/, a,z;G) is computed in the Parry's 
limit as 

4(3 ~" )/v/2 n - r 1 E r (^^g+p) c .^" 



r(|) Al!Ko 



(19) 



4(3-^/2^1/2 ^ 



.4 



1 

2^ 



»+oo 



alk 



r(¥^ 2 + ^ 2 )) 

(G 2 + P) (3-l)/2 



The second method to derive the Ewald method for quasi-two dimensional systems 
follows exactly the same derivation as the one done in section 2 for tridimensional 
systems. We begin by applying Eq. (jSJ to the lattice sums ( TT21) and we split the integral 
into two contributions introducing a. Then, we have 

v2 



ir+S, 



r(I) 



dt 



-Z 2 t 



E 



3 -|s+5 n | 2 t 



The Poisson-Jacobi identity in two dimensions, with IZ2 the reciprocal lattice, 



-|s+5nl 2 * = I ^ 



]T e^exp 



S n G&l 2 
is applied to the second contribution, thus we obtain 



G]_l_ 



(20) 



(21) 



\(f) =^$g ) (r / ,a;|r + 5 n |) 



1 ^ Q S 



E 



Sn 

1 7T 



2 + _G 



(22) 



r(2M./ t M) 
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Power 



Contributions depending 
on z 



Divergences 



T) > 2 

j] = 4 

V = 3 

77 = 2 

v = i 



|(^) 7 (|-l,aV 



r(f)4 



2tt 1 



erf (a | z |) 



A 1 
2tt 



| z I erf(a | z |) + 



None 



None 



None 



aV)] 


27T 

— - me 
A 


2 2 ~ 

e -a z 


2^1 




.4 e 



Table 2. Non trivial contributions depending on z in quasi- two dimensional 
Ewald summations and asymptotic behaviour for G = 0. 



The last contribution in Eq.( l22l is the contribution for G — 0, it corresponds to the 
contribution Tq_^ in Eq.(fl~6l). With this derivation, we may write (f> v {r) as in Eq.(TT2l) 



with 



and 



(„) _^7T I- zH 

j g=o- r(p ' -"-^ e 





2 ' 


1 


7T 


r(f 


M 


i 


7T 



(24) 



r(f) a v 4 ; y GV4a2 



Few analytical forms of Tq_ q , with their dependence on z and their asymptotic 
expension in e (if n < 2) are given in Table [2] for few integer values of r\. 
Two particular cases can easily be obtained with Eq. (I24j) : z — 0, that corresponds to 
the Ewald method for a two dimensional system (see also [13]) ; thus we have 

*«?(„,a;G) = •g»(,,a,0;G) = ~ (f )'""' r (l - §, (25) 
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and the limit a — > oo, that gives the Nijboer-de Wette representation for inverse power 
law interactions in quasi-two dimensional systems [Ml 05] and then we found 

*25, (i7,s;G) = <S> ( § 2) ( v ,oc,z ] G) 

| * K (| _ x) (G | * |) 



r(2) a V 2 



The general case is more complicated to obtain, it is necessary to use incomplete Bessel 
or generalized incomplete gamma functions [32| |34"1 135] . Following the definition given 
by Harris in ref . [32] , the incomplete Bessel function is defined as 



oo 



K*(*,y) = / (27) 
Then, with the second integral in Eqs. (l24[) . we find 



9%"(r,, a ,z;G) = WT K (| _ 1} ( — ,«V ) (2* 

-2> 



According to the connection between the generalized incomplete gamma function and 
the incomplete Bessel function, defined as (Eq.(6) of ref. [32] ) 

K v (x,y) = x v T (-v,x\xy) (29) 

Eq.( l28j) can also be written as 

If rj — 1, the inverse power law interaction corresponds to the Coulomb interaction ; one 
can easily verify that Eqs. (l28|) and (1301) allow to recover the Ewald sums for Coulomb 
interactions in quasi-two dimensional systems [231 133 SOI EJ E2] (see also Theorem 7 
in ref. [33]). 

On Table El we give some analytical formulas for the computation of the reciprocal 
space contributions in quasi-two dimensional Ewald summations for inverse power 
law interactions. For numerical implementations, convenient algorithms for the 
computations of the incomplete Bessel functions are needed (see for instance refs. [36l 

EI). 

To close this section, we want to outline that the equivalence between both derivations 
allows to obtain a new analytical relation for the Fourier transform of incomplete gamma 
function. From Eq. (Tl9l) . we find 

1 r d Ji^MG 2 + ^)) ( 2tt )&HP / G 2 22 



1 _ v ( 1 _r L G\G 2 z 2 



(31) 



2 V ^FG( 2 -") V 2 '4a 2 ' 4 
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Power 
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Eq.(20) in ref.f^2l 
(Qt?i( 2 ) 1S a Legendre function (cf. section 8.7 of ref. |38|) 
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2 _ y2 J J 


T) = 3 


* 1 r -a, t ( G \ a* . ( G 

e eric 1 — olz 1 — e eric I - 

A | 2 | L V 2a / V 2a 


- cez^ 




Eq.(28) in ref.f34l 
and also in Ewald methods for dipolar interactions. 


7, = 2 


: [ko(gud e -^+xh £ (^ 2 -y 2 ) m q 

A [ _~ Q m! 
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Eq.(lD) in ref.f32l 


7, = 1 
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e eric 1 — a2 1 + e crfc 1 + 

A G L V 2a / V 2a 


otzj J 




Theorem 7 in ref. I34| 
and Ewald methods for Coulomb interactions. 



Table 3. Analytical formulas for the reciprocal space contributions in quasi- 
two dimensional Ewald summations for inverse power law interactions l/r r ' 
for some values of r\. The column 'References' gives some sources useful for 
analytical computations of ^q 2 \w, a, z;G). As shown with Eq. (|12p . the real 
space contribution to Ewald summations is given in Table [TJ 



If one takes a = 1 and G = 0, we recover the Fourier transform of the incomplete 
gamma function for v > 1/2, as 



where we have used the definition of the incomplete gamma function 7(/x, y) 

7 (/i, y) = I" dx x^-V e~ x = y" KJO, y) (33) 



Conversely, an analytical demonstration of Eq. (|3T|) based on the properties of incomplete 
and generalized incomplete gamma functions will demonstrate the equivalence of both 
derivations for inverse power law interactions (already shown for Coulomb, Yukawa and 
Dipolar interactions [241 l25l H8]). Such a pure analytical derivation of Eq. (j3TT) has not 
yet been achieved [33] . 

The equivalence between both derivations may also be justified directly by the relation 
between Poisson-Jacobi identities in three and two dimensions. For one dimensional 
periodicity, the Poisson-Jacobi is [15] 
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then, multiplying Eq.( l2ip with Eq. ()34p . we obtain the Poisson-Jacobi for three 
dimensions Eq.flS]). This agrees with the separations done in Eqs. ()13|14j) . 

4. Electroneutrality and IR-divergences. 

The IR-divergences, that occur in Ewald sums of inverse power-law interactions when 
n < d, with d the dimension of the periodicity, are related to the conditional convergence 
properties of the lattice sums Eq.flTJ). These divergences are an artefact due to the 
approximation of the finite lattice sum by an infinite lattice sum (Ewald sum) [121 EH] ■ 
For coulomb interactions these divergences in the infinite lattice sums are cancelled in 
the computation of energy and forces when the system in the basic cell is electroneutral. 
Two classical models are used with coulomb interactions : the restricted primitive model 
of electrolytes (RPM) and the one component plasma model (OCP) [471 SHI US]- The 
RPM is made of 2 No hard spheres that carry electric point charge at the centre of 
the sphere, N particles carry a charge +Q and the other N particles carry a charge 
—Q. We define the 77-RPM similarly to the standard RPM, but with interaction between 
charges as inverse power law interactions. In the OCP model, there are N point particles 
carrying all the same charge Q and a constant volume density of charge po- in the rj- 
OCP model, the system is defined as the standard OCP model and interactions between 
point particles are inverse power law interactions. 

In the following, we first compute the energy of a system of N point particles carrying 
pseudo-charges Qi in a system with periodic boundary conditions in all three directions 
of the space. Then, in subsection 4.1, we make the same computation but for a quasi- two 
dimensional system where electroneutrality is fulfilled with the help of a constant planar 
surface density of pseudo-charge. In subsection 4.2, we apply the analytical results of 
section 3. to a monolayer system made of N point particles carrying the same pseudo- 
charge Q, this model is a generalization of the two dimensional One-Component plasma 
(77-OCP). Few computations using Monte-Carlo simulations for monolayers of 77- RPM 
and 77-OCP systems are reported in subsection 4.2. 
Let the interaction energy between two pseudo-charges be 



where ro is a typical length defined by the geometry of particles, for instance the diameter 
of hard spheres in the 77- RPM model. In the following, all length are measured in unit of 
r (i.e. we set r = 1). The energy of the system with tridimensional periodic boundary 
conditions is given by 




(35) 




(36) 



i=l j^i 



i=l 
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with the contribution of the interaction of a charge with its own periodic images. 
This quantity is defined and computed as for <f) v (r), thus we find 

1 E^ 3) (^«;i^ni) + E $ i 3) (^ a ; fc ) 

(37) 



^ = £ 



-a 71 



1 7T 3 / 2 r* alt 



r/r(f) r(a) v y ^- v)/2 

The last contribution in the previous equation has the same IR-divergence than in 
Eq.()9]). When it is grouped with the other contributions of Eq. fl36|) . we have 




t(5-,)/ 2 - ( 38 ) 

that is cancelled if Ylf=i Qi = 0- ^ would be more rigorous to consider that this equation 
is an indeterminate of the form x oo ; it corresponds to the property of conditional 
convergence of lattice sums for inverse power-law interactions, this is related to the 
macroscopic boundary condition and to the particularity of the systems (see for instance 
ref. |46[ H~2] for Coulomb interaction, rj = 1). 
Then the energy can be written as 

^ N N 



2 

j=i i=i j jn 



(39) 



N 



E^ ex P {jk.ri) 



N 



»r(f) 



where the prime in the first contribution indicates that the term i = j is not included 
when L n = Lq = 0. We have also made a factoring of the reciprocal space contributions 
as 1-particle summations, this is an important property for the numerical efficiency of the 
summations as it allows to the computational time needed for obtaining the reciprocal 
contribution to the energy to scale as iV and not as n 2 . One may verify easily that for 
n — 1, we recover the Ewald summations for Coulomb interactions, Eq. fl39|) allows to 
obtain the energy of the 77-RPM tridimensional systems. 

In the 7/-One Component Plasma model (77-OCP) all point-particles carry the same 
charge Q and electroneutrality is achieved by the constant volume charge density po 
such as 

NQ + p V = (40) 

For this system, one may show easily that the IR-divergence is cancelled as in Eq. flHS]) . 
written as 

1 7r 3/2 r a2 dt 

(NQ + PoVf 77^ = (41) 



r(a) v v ^ KU 1 J t( 5 -")/ 2 
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The energy E{rj-OGP) is given by using Eq. ()39p with Qi = Q, for < n < 3, we obtain 



7r 3 /V"- 3 ) N 2 Q 2 



E ^- 0CP) = E «-(3^WW)^ (42) 

and for 77 = 3, we have 

N 2 2 

£(3-0CP) = E 3 - 2vr— ^- In a 2 (43) 

In Eqs. fl42|43p . the contributions that scale as N 2 stem from the interaction of particles 
with the neutralizing background. For n > 3, no IR-divergence is present and the energy 
of the system can be obtained from Eq.( l39i) . 

4-1. Electroneutrality in quasi-two dimensional systems. 

To extend the tridimensional 77-OCP model to quasi-two dimensional systems, one has 
to choose a particular neutralizing background to fulfill electroneutrality ; a constant 
volume charge density po is inconsistent with the symmetry of quasi-two dimensional 
systems. Several choices can be done [25] ; in the following, we choose the neutralizing 
background as a constant planar surface density a of pseudo- charges. With this choice, 
the electroneutrality of the system reads as 

NQ + a A = (44) 

and the energy of the 77-OCP, for < n < 2, can be computed as 

Y2 



E(V-0CP) =g-EEE / . r .. + 1 5w |n +^oE/, ^ |r _ s 1 +Snh 
1 i= i j= i Sn I r » + *n I . =1 JS Sn \r t s + bn\ 

(45) 

+ °4 I ds [ ds'Y— 1 —— 

2j So Js £|-'-« + S„|i 



5... 

To compute the lattice sums, we apply the methods of section 2 and 3 ; with the 
notations of Eq. ffTBl . for < n < 2, we find 

/02 N N ( V (v-2) \ 

£( "" 0CP, " T £ £ £ <* - 1 - + S " " + 3 (Fi^) ^ 

4 E E E «. %; g) - 2^7^% £ 1 * i' 2 -"* («) 
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and for 77 = 2, we have 



£(2-ocp) = ^EEE ( 2 ' «; 1 r « + 5 ^ I) 



AT JV 



TriVQ 2 / nl a 2 A 
' C + 2 In a 



i=i i=i 5 



2A 



7T 



77, 



GMo i=1 J =1 



^ Q2) (2,«,^;G)- 



nCN 2 Q 2 
2A 



(47) 



2A 

i=i j^t 



Ei (« 2 4) 



Win 



E ln * 



1=1 



with C the Euler's constant and Ei(x) the exponential integral. 

It is worthwhile to note that for quasi-two dimensional systems the factoring of the 
reciprocal part of the energy into 1-particle summations cannot be achieved because of 
the complicated dependence on z in , a, z; G). 



4-2. Monolayers and few numerical results. 

When all particles are confined in a plane, the system is a monolayer and one can factor 
the reciprocal part of the energy into 1-particle summations. For a system of N point 
particles carrying a pseudo-charge Qi, satisfying to the electroneutrality, one finds easily 



N N N 

E t ] = o E E QiQi E '** fo. «; 1 + 5 * i) - E Q 



i=i j=i 



(!) U 



N 



exp (jG.Si 



i=i 



For the ?7-One Component Plasma model with all point particles confined in the plan 
with the constant surface charge density, the energy is given by 



E,, 



0? 

2 



N N 



EEE'^ tv,<x;\s ij + S n \) 



a- 



i=i j=i s 



n 



»?r(l) 



NQ 2 



(49) 



+q 2 X>gW;G ! ) 

G^o 



N 



1=1 



7T a 



(>7-2) 



2/o2 



A(2-^)r(2) 



iv^g 



As outlined in the introduction, to recover the correct Ewald method for coulomb 
interaction in two dimensions, one may not take the limit n — > in Eqs. (148f49j) . but 
one must use Eq.([3]). For coulomb interactions (77 = 1), monolayers of the restricted 
primitive model of electrolytes |49] and of the one component plasma model [48] have 
been studied previously. 
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77-RPM 






77-OCP 


n 


I3U/N 




V 


/3U/N 


0.25 


-10.47(3) 


5 x 10 2 


0.5 


-10.76(2) 


0.5 


-11.36(3) 


5 x 10 


1.0+ 


-16.67(3) 


0.75 


-12.29(3) 


5 x 10 2 


1.5 


-33.04(3) 


1.0* 


-13.23(3) 








1.25 


-14.21(5) 


5 x 10 3 


0.5 


-35.26(3) 


1.5 


-15.2(2) 


5 x 10 3 


1.0* 


-54.18(3) 


1.75 


-16.3(3) 


5 x 10 3 


1.5 


-105.97(3) 



t In this table, e is defined as in ref.[35] : e = 2tt pQ A /T 2 . 
* cf. ref.gS], Table n ; Line 15 - 
* cf. ref.gg], Table II. 

Table 4. Average energies for 77-RPM and 77-OCP models in Monte Carlo 
computations. The numbers in brackets give the accuracy on the last digit of 
the averages. For the 77-RPM, the average energy is defined by U =< E^ >mC 
with given by Eq, (|48p ; for the computations reported in this table, one 

has : N = 2N = 1024, p = 2N /A = 0.6, Q 2 = 20 (it corresponds to /3* = 20 
in notations of ref.[39]). For the ?7-OCP model, the average energy is defined 
as U =< E v >mc with Eq.flU]) ; one has : N = 1024, irp = ttN/A = 1.0, 
Q ~ 3.98 for e = 5 x 10 2 and Q ~ 7.071 for e = 5 x 10 3 . 

In the following, we report some preliminaries numerical results obtained by Monte- 
Carlo simulations in the canonical ensemble (NAT) [50]. For 77-RPM, the computations 
are done in a square box with a fixed shape, while for 77-OCP monolayer models the 
surface A of the simulation box is fixed, but the shape of the basic cell is allowed 
to fluctuate [51]. Periodic boundary conditions and the minimum image convention 
are applied. For both models, average energies are computed from Eqs.(j4HH9j), pair 
correlations functions are also exactly computed with the same definitions and methods 
as in refs. [48l l50| 1491 [5ll [52]. Voronoi constructions and cells for the 77-OCP monolayers 
are exactly computed as in ref . [52"] . 

In Table HI we report some average energies computed with Monte-Carlo simulations for 
the r/-RPM and 7/-OCP models. The values obtained with 77 = 1.0 do well agree with 
previous results obtained in refs. [HH 08] . 

On Figured], we show some snapshots for 77-RPM monolayers (a-b) and 77-OCP models 
(c-d). Snapshots with 77 = 1.0 agree with the structure found in systems with coulomb 
interactions [H1H9], while for 77 > 1.0 short ranged order is more marked. This finding 
is enforced with the shape of pair correlation functions given in Figure [2] 
On Figure [2] (a), pair correlation functions in the 77-RPM monolayers between hard 
spheres that carry positive charge are represented. For 77 = 1.0, g++(r) is exactly the 
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Figure 1. Snapshots of 7/-RPM and n-OGP monolayer models. The sides of 
the simulation box are represented by thick black lines, periodic boundary 
conditions are applied, (a-b) : 7/-RPM, blue disks are particles carrying a 
negative charge (— Q) and red disks a positive charge (+Q). For both snapshots, 
N = 2Nq = 1024, p = 2N /A = 0.6, Q 2 = 20 (cf. Table j ; (a) r) = 1.0 and (b) 
7] = 1.75. (c-d) : ?/-OCP model with Voronoi construction, for both snapshots 
: N = 1024, up = ttN/A = 1.0 and e = 5 x 10 3 (cf. Table [J). Voronoi cells 
with four sides are represented in yellow, those with five sides are represented 
in green, those with six sides in white, those with seven sides in red and those 
with eight sides in blue. Point particles are represented by a small dot. (c) 
7] = 1.0 and (d) n = 1.5. 

same as the one obtained for the standard RPM model reported on FIG. 10 (a) of ref. [32]. 
As already shown on the snapshot of Figure (TJ when rj > 1 the short ranged order and 
correlation between particles are more pronounced than with coulomb interaction ; they 
are less if rj < 1. A similar behavior is found in ?7-OCP monolayers ; on Figure W[b), 
we report pair correlation functions between point particles in 77-OCP monolayers. For 
rj = 1, g(r) is exactly the same as in FIG. 2 of ref. [48] . 
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2.5 I — 1 — ' — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — ' — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — ' — 1 — 1 — i — 1 — ' — 1 — i — 1 — 1 — 1 — i — 1 — 1 — r 




Figure 2. Pair correlation functions for ?]-RPM and ?7-0CP models. The 
thermodynamic states of each system are the same as those reported on Table 
[U For each system, the values of n are given in legends, (a) Pair correlation 
functions g++(r) between positive ions in the 77-RPM model ; for 77 = 1.0, 
g++(r) can be compared with the standard RPM model (see FIG. 10(a) of 
ref.[l9]). (b) Pair correlation functions g(r) between point particles in the n- 
OCP model with e = 2irpQ 4 /T 2 = 5 x 10 3 ; for n = 1.0, g(r) is the same as the 
one given in the FIG. 2 of ref . |48| . 
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5. Discussion. 

From the potential <p v (f) computed with the Ewald method for tridimensional, quasi- 
two dimensional or two dimensional systems derived in section 2 and 3, one may easily 
obtain any physical quantity related to the potential in a similar way as it is done for 
the energy. 

In ref. [31], Johnson and Ranganathan have proposed a generalized approach to Ewald 
sums for diverse long ranged central potentials, including inverse power law potential as 
(f)(r) = r - ( 1+<5 ) with < 5 < 1. Their derivation is based on the split of the interaction 
potential as (cf.Eq.(6) of ref . [31] ) 

0( r ) = (1 - /(r))0(r) + f{r)<p{r) = (1 - erf^)) <j>(r) + erf(ar M )</>(r) (50) 

where a and /x are chosen conveniently. In lattice sums, the first contribution in the 
right hand side of Eq. (1501) is evaluated as the real space contribution and the second as 
the reciprocal space contribution by using the Fourier transform of f(r). 
In the present paper, the derivation of Ewald methods for inverse power law interactions 
is done by using the Poisson- Jacobi identities ; according to the results for tridimensional 
systems in Eqs. (l9|lip . it corresponds to a choice of the screening function as f(r) = 
7(77/2, a 2 r 2 )/T(rj/2). For coulomb potential (77 = 1, or 5 = in notations of ref. [31] ) 
both choices are strictly equivalent, for other values of n formulas differ since the choice 
of the screening of pseudo-charges differ ; however, both derivations should lead to the 
same numerical results [261 [29] . 

The preliminary numerical results for 77-RPM and 77-OCP monolayers, given in section 
4., show that the Coulomb potential case is correctly reproduced by inverse power law 
interactions with 77 = 1 when compared to previous results [4*8l 149] . A longer numerical 
study of these systems is ongoing. 
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